Non-linear mechanical response of the Red Blood Cell 
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We measure the dynamical mechanical properties of human red blood cells. Single cell response is 
measured with optical tweezers. We investigate both the stress relaxation following a fast deforma- 
tion, and the effect of varying the strain rate. We find a power law decay of the stress as a function 
of time, down to a plateau stress, and a power law increase of the cell's elasticity as a function 
of the strain rate. Interestingly, the exponents of these quantities violate the linear superposition 
principle, indicating a nonlinear response. We propose that this is due to breaking of a fraction of 
the crosslinks during the deformation process. The Soft Glassy Rheology Model accounts for the 
relation between the exponents we observe experimentally. This picture is consistent with recent 
models of bond remodeling in the red blood cell's molecular structure. Our results imply that the 
blood cell's mechanical behavior depends critically on the deformation process. 

PACS numbers: 87.68.+Z, 83.60.-a,87.16.-b,87.17.-d,87.80.Cc 



The human red blood cell is a biological structure of 
relative simplicity: it lacks a nucleus and intra membrane 
organelles. Its components are well characterized 
and they can be summarized as a bilayer membrane, cou- 
pled to a thin tenuous cytoskeleton of spectrin filaments 
via two complexes of a few different proteins (ankyrin, 
band 3 and band 4.1). This outer membrane, held un- 
der tension by the cortical cytoskeleton, encloses a so- 
lution of dense haemoglobin. Despite its simplicity, the 
properties of this structure have puzzled researchers for 
decades. Investigations have focussed on two, related as- 
pects: explaining the biconcave discocyte shape that is 
found under physiological conditions (and a whole array 
of other morphologies that are found under perturbation 
or disease); measuring the mechanical properties of the 
structure. Regarding the shape, significant progress has 
been made recently showing how, by carefully balancing 
the membrane bending elasticity with a tension from the 
underlying spectrin scaffold, the discocyte shape emerges 
as the equilibrium solution, and the cup-shaped stoma- 
tocyte is the result of a perturbation that decreases the 
membrane area relative to the underlying skeleton[3| . On 
the issue of mechanics and dynamics of the structure, our 
understanding is still very limited, despite the great im- 
portance of these parameters. Large deformations of the 
cell are involved in blood flow through thin capillaries, 
and the cells in suspension more generally determine the 
peripheral blood rheology in healthy and diseased condi- 
tions [4] . The red blood cell has variously been described 
as either liquid or soHd 

In this Letter we investigate the dynamical mechanical 
properties of the red blood cell by measuring its vis- 
coelasticity over a range of environmental conditions. We 
focus on the effect of the cell age, of the deformation 
dynamics, and the correlation between the cell modulus 
and its shape. We relate the findings to phenomeno- 
logical models for flow in kinetically arrested systems, 
and to recent models describing the metabolic activity 
of the RBC, where the consumption of ATP controls 
the stiffness of the elastic network d, 0, Hi- The nor- 



mal biconcave human RBC has approximately 25% more 
membrane surface area than the minimum required to 
enclose its volume, allowing deformations without exten- 
sion of the membrane [9j. Therefore the mechanical re- 
sponse of the cell is related only to its bending and shear 
elastic moduli, and its membrane and bulk viscosities. 
These quantities have been measured with various tech- 
niques, in particular micropipette aspiration, shape flick- 
ering, deformation in flow, and optical tweezers. Existing 
measurements give very different values. For example, 
independent experiments using micropipette aspiration 
yield a shear modulus of between 6 and lO^N/m [lo} . 
Optical tweezers measurements report a wide range of 
shear modulus, but the discrepancies between these re- 
sults are in part the result of applying different geomet- 
rical models to extract the modulus from the measured 
forces: 2.5±0.4/iN/m assuming a flat disk geome try [Tlj : 
200/iN/m assuming a spherical initial shape [H, [HI; 11 
to 18/iN/m comparing to a flnite element simulation of 
the deformation p3|- A satisfactory approximation for 
the deformation geometry is still not available, and is not 
the topic of this Letter. We choose instead to report the 
directly measured mechanical properties, showing that 
even the model-independent dynamical stiffness of the 
cell is strongly dependent on the deformation protocol. 
Our key result is that the RBCs response cannot be ex- 
plained within the framework of linear viscoelasticity. 
Fresh blood was drawn from a healthy volunteer donor, 
and diluted in phosphate-buffered saline (PBS) with 
acid citrate dextrose (Sigma C3821) and Img/ml bovine 
serum albumin (BSA) (Sigma A4503) at pII7.4. For 
studying the effect of aging time after drawing, cells were 
incubated with glucose-free PBS either in a hot bath 
plate at 37°C for 24h with 10~^ Penicillin-Streptomycin 
solution (GIBCO, Invitrogen 15070063) to prevent micro- 
bial growth, or at 4°C for 24h. A flrst control group was 
formed by adding ImM MgATP (Sigma A9189), after 
the 24h incubation, to the the RBCs in glucose-free PBS 
suspension. A second control group was incubated with 
lOOmg/dl glucose in PBS buffer. These protocols provide 
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samples that are a mixture of discocyte and stomatocyte 
cells, with an increasing fraction of discocytes after 1 
day's incubation. 

The optical tweezers setup consists of a laser (IPG 
Photonics, PYL-1-1064-LP, A=1064nm, P^aa;=l.lW) fo- 
cused through a water immersion objective (Zeiss, Achro- 
plan IR 63x/0.90 W) trapping from below. The laser 
beam is steered via a pair of acousto-optic deflectors (AA 
Opto-Electronic, AA.DTS.XY-250@1064nm) controlled 
by custom built electronics allowing multiple trap gener- 
ation with sub nanometer position resolution. The trap- 
ping potential is locally described by a harmonic spring, 
and the trap stiffness was calibrated by measuring the 
thermal displacements of a trapped bead. The trapping 
stiffness at maximum power when trapping two beads (as 
in this work) is fc|;j.g^p=44pN//im on each bead. The sam- 
ple is illuminated with a halogen lamp and is observed 
in bright field with a fast CMOS camera (Allied Vision 
Technologies, Marhn F-131B). 

Carboxylated silica beads of 5.0/im diameter (Bangs 
Labs) were washed in Mes buffer (Sigma M1317). They 
were then functionalized by resuspcnding in Mes buffer 
and incubating for 8h at 37°C with Lectin (0.25mg/ml) 
(Sigma L9640) and EDC (Sigma E1769) (4mg/ml). This 
made them very sticky to the RBC, probably through 
binding to glicoproteins or sugars on the outer side of 
the RBC membrane. To prevent the RBCs from sticking 
on the glass surfaces of the chamber, BSA was coated 
on the slide glass. Using the tweezers, two beads are 
brought to an RBC and attached diametrically across a 
cell, which is then floated well above the glass slide sur- 
face (at about 10 times the bead diameter) to minimize 
any hydrodynamic drag from the solid surface. The cells 
are quite monodisperse in size, with an initial cell length 
(the diameter) Lq ~ 8/.tm. The area of contact between 
the bead and the cell varies in the range 3-4. 5^m^, and 
we find no correlation of the patch area with any of the 
experimental results. This arrangement is drawn in Fig- 




FIG. 1: Diagram of experimental setup. A cross marks the 
positions of the laser trap and a star marks the center of the 
bead, obtained via image analysis, (a)-(d) Decays of the mea- 
sured time-dependent cell stiffness K{t) as a function of time. 
The plateau value has been subtracted. This data shows a 
power law decay with exponents ~3/4. 



ure [TJ During all the deformation protocols, one laser 
trap is kept fixed at its initial position, while the other 
trap is moved away by a distance Alaser. We focus on 
tracking the bead in the immobile trap, measuring the 
difference Ax between the bead and the known laser po- 
sition using image analysis code written in Matlab. From 
the displacement Ax, the stretching force is calculated. 
The cell is elongated by a distance AL = Alaser — 2Ax, 
and we define the strain as 7 = AL/Lq. The resolu- 
tion of bead position via image analysis on each frame is 
around 5nm, which translates into ±0.22pN of force res- 
olution. This is significantly less than the random fluc- 
tuations caused by thermal noise. The cell stretching is 
recorded at ~60 frames per second, and having checked 
that we could not observe a variation between measure- 
ments done at room temperature and at 37°C, we report 
on measurements made at room temperature. Tempera- 
ture is known from previous studies to have only a slight 
influence on the mechanics of red blood cells [3, EB] ■ 
We perform stress relaxation experiments by moving just 
one trap (i.e. moving one bead) by Alaser=3.5/^m at 
20/im/s, and monitoring the force acting on the station- 
ary bead. The force is observed to decay towards a 
plateau value, which is reached within less than half a 
minute. We define the time-dependent stiffness of the 
red blood cell to be K{t) = F{t)/AL{t). The whole set 
of experiments is fit very well by the 3-parameter power 
law function: 

K{t) = Koo + AK(i) = Koo + AKo(iAo)"", (1) 

where we fix to = Is to obtain a dimensionless time. We 
use the first 20 seconds of the decay to fit these values, 
where there is less noise in the data, but it can be seen 
that the power law decay holds for longer times as well. 
In contrast, a single exponential fails to fit the data and 
even a stretched exponential (4-parameters) gives a poor 
fit. Figure [1] shows examples of the time relaxation of 
AK{t) plotted separately for each of discocyte and cup 
shaped stomatocyte cells, fresh and after one day of in- 
cubation. In all cases the stiffness is fit over two decades 
in time by the form of Eq. [T] The power law exponents 
are the same within experimental error, the mean of all 
the power law exponents measured is a=0.75±0.16 for 
fresh cells and a=0.82±0.09 for 1-day old cells. A to- 
tal of 11 different cells were measured this way. The 
value of the ratio AFq/Foo, i.e. the fraction of the stress 
that decays over time, also appears independent of shape 
and cell age, and is between 0.55 and 0.75. The fact 
that the data is described so well by Eq. [1] implies that 
there are two components to the cell elastic modulus: 
a time-dependent and a long-time equilibrium modulus. 
Power-law relaxations are well known in systems near the 
gelation point and in entangled polymer solutions [13], 
however we anticipate here that we will discuss how for 
the red blood cell the form of Eq. [1] does not originate 
from polymer relaxation nor gel cluster dynamics. To 
investigate further the strain rate dependence of the cell 
modulus, we perform a triangle wave deformation exper- 
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iment using speeds of 20, 5, 1 and 0.2^m/s, and 3.5^ni 
laser trap displacement Figure [5{a) shows that the 
force measured over time is approximately proportional 
to the applied strain. The slight mismatch between the 
applied strain and the stress can be seen more clearly in 
Figure El^b) as a hysteresis loop. There is a slight dif- 
ference between the onset of the first cycle and all the 
following ones: From observing many instances of this 
protocol we believe that this is due to the cell rearrang- 
ing its position precisely along the axis of the strain ex- 
periment. It is clear from this figure that there is no 
further evolution of the stress strain curves over the run, 
which lasts typically tens of cycles. The cell's response is 
predominantly elastic, and the gradient of the force vs. 
elongation data as in Figure [2llb) , taken over the exten- 
sion, gives the stiffness k of the cell. This stiffness is an 
effective property of the cell, a valid modulus for describ- 
ing the cell's response in physiological processes such as 
squeezing through thin capillaries. Our results show that 
there can be a large (factor of three) difference between 
the dynamic and static stiffness, as measured for example 
by micropipette. The experimental speed of deformation 
can be converted into an approximate strain rate 7 by 
dividing the bead velocity by the initial cell length (Lq). 
The cell stiffness is plotted in Figure [3] against the strain 
rate. These quantities are seen to be approximately re- 
lated through the function: 

k(7) - ^0 7^ (2) 

with /3 ~ J. 

The hysteresis in the stress strain curves is a measure of 
energy dissipation in the cell, and is shown in Figure [31 
In order to compare results from different strain ampli- 
tudes, the dissipated energy is presented as the work per 
cycle, normalized by the square of the strain amplitude 
Imax- The dissipated energy follows a power law as a 
function of strain rate, with an exponent that we assume 
is the same /? as in Eq. [21 In all the cells the strain rate 
dependence is well below the linear relation expected for 
a purely viscous system. The exponents fitted here are 
f3 =0.25 for fresh discocytes and (3 =0.33 for discocytes 
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FIG. 2: (a) The dashed line is the elongation AL of the 
red blood cell subject to a triangle wave trap displacement, 
and the solid line is Ax of the bead in the stationary trap, 
(b) Force vs. elongation of the cell AL, for the same data 
as in panel (a). Solid line is a linear fit to the force during 
extension, the area inside the hysteresis loop is the energy 
dissipated per cycle. 



after one day's incubation. These values have a higher 
precision than the stress relaxation exponents because 
the triangle wave cycle is repeated many times. (26| In 
linear viscoelasticity, the dependence of the modulus on 
the strain rate is related to the form of the stress relax- 
ation via a modified Laplace transform [iTj . In the case 
of a power law decay with K{t) — k^o we expect 

k(7) — ko ~ 7". This is very clearly not the case here, 
the exponent (3 is much smaller than a. This points to 
the fact that even though our strain amplitudes are not 
very large, only of the order of 20%, our experiment is 
in a nonlinear regime. This finding can be understood as 
follows: the power law relaxation is the manifestation of 
a system where at least a fraction of the bonds can break 
under stress and quickly reform, essentially remodeling 
the skeleton's network. There is an analogy between this 
molecular-scale picture and the more general situation of 
the rheology of a system where the components have to 
overcome potential barriers in order to flow. This condi- 
tion is addressed by the soft glassy rheology (SGR) model 
which, assuming a distribution of energy wells of different 
depths, provides a framework for calculating the nonlin- 
ear response of the system [l8|. It contains one principal 
parameter, x, which is the ratio between the available 
energy and the mean well depth. For x < 1 the model is 
in a glassy phase. The shear modulus is predicted to de- 
cay with time as C? and the stress under constant 
strain rate \s a ~ Cy ^ 7^^^. The red blood cell power 
laws and exponents appears to be well described by the 
SGR model with a value of a; ~ 3/4. In the red blood 
cell it is likely that the potential barriers to be overcome 
are the energies for releasing a spectrin filament from a 
crosslink. The filament would then re-bond in a con- 
figuration of lower stress. The presence of a long lived 
residual stress in the stress relaxation implies either the 
presence of a fraction of permanent bonds, or of a ki- 
netically arrested state with residual stress, leading to a 
deviation from a purely power law scaling at low strain 
rates which is visible in the stiffness data of Figure [31 
In the absence of external stresses it is known that ATP 
is required to open the spectrin-actin linker, whose bond 
energy is around Tk^T, i.e. 3xI0-20j [11, HO]. Based 
on the known mesh size of the cytoskeleton (between 80- 
lOOnm) and the cell surface area, we can estimate the 
order of magnitude of the number of bonds to be 10^. 
Given that the dissipated energy per cycle is of the order 
of 10^^^ J this implies that either only a small fraction of 
the bonds, around 3%, is broken during a cycle, or that 
the bond energies under stress are reduced considerably, 
or that ATP and not solely mechanical energy is being 
used to break the bonds. Of these possibilities the agree- 
ment with the soft glassy rheology observed here suggests 
that remodeling of the linkers under stress does not in- 
volve the full cost of breaking the spectrin-actin linker. 
While the features of the dynamical response of all the 
cells are the same, we observe differences in the magni- 
tude of the response depending on the shape of cell and 
the incubation protocol. The strongest influence, lead- 
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FIG. 3: Strain rate dependence of the RBC stiffness k and 
the energy dissipation. Markers indicate the different condi- 
tions studied here: (t) fresh-Discotic; (v) fresh-Stomatocyte; 
(•) f day-Discotic; (o) Iday-Stomatocyte; (■) f day-discotic 
+MgATP. 

ing to a decrease of the cell modulus of almost a factor 
of 2 after 1 day of incubation, is the time from when 
the cells are drawn. Secondary to this aging, we ob- 
serve a correlation between the stiffness of the cell and 
its shape. The cup-shaped stomatocytes are almost twice 
as stiff as the discotic cells. Models proposed recently de- 
scribe the active remodelling of the RBC's cytoskeletal 
network [1, [2^ . ATP is expected to be strongly related 
to this activity, at least under absence of an external 
force ft is of interest that flickering of the membrane 
at low frequencies was shown to depend on the intracel- 



lular MgATP [2l[ , and the flickering has also been shown 
to decrease with cell aging [l^. Through the 1-day in- 
cubation protocol the ATP concentration inside the cell 
is depleted, and according to these references this would 
lead to stiffer cells. Our results showing softening seem 
to imply an opposite effect, one possibility is that under 
reduction of ATP the cell loses some of its ability to re- 
form the bonds that break under deformation. We find 
no appreciable difference between aging in the absence of 
glucose and in physiological glucose concentrations (data 
not shown). We find a strong effect adding MgATP to an 
aged cell: MgATP causes an increase of both the stiffness 
and dissipation for the discotic cells, bringing the values 
close to those of freshly drawn cells. 
The rate dependence of the cell modulus measured in this 
work can be also be compared to recent data obtained by 
monitoring the deformation induced via magnetic twist- 
ing cytometry (MTC). This is essentially a creep exper- 
iment performed on a pivoting rnagnetic bead bound to 
the outside of a red blood cell [22|. In MTC there is 
almost certainly a large influence of the membrane bend- 
ing elasticity, and the data in ref. [1^ is dominated by 
a source of elasticity that is not found in our measure- 
ments with optical tweezers. Our data confirms that the 
dynamics in the red blood cell, as also pointed out by 
Suresh and collaborators, is intrinsically free of a char- 
acteristic timescale [225 . In addition we have shown here 
that the dynamics is not linear, and therefore depends on 
the deformation protocol, even for small deformations. 
This may be relevant to other cell types. Our findings 
here are in contrast to many previous investigations of 
the red blood cell [ll|, [H, [2^ , where limited datasets led 
to simplistic conclusions regarding the dynamics. 
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